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BRUSHLESS MACHINE CONTROL 

CROSS-REFERENCE TO RELATED APPLICATION 

5 The subject matter of this application is related to the subject matter of Application 
No. GB 9929994.3, filed December 17, 1999, priority to which is claimed under 35 
U.S.C. § 1 19 and which is incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

10 

1. Field of the Invention 

This invention relates to the control of electronically switched, brushless machines, 
such as switched reluctance machines, permanent magnet machines and hybrids 
1 5 thereof. 

2. Description of Related Art 

The typical switched reluctance machine, for example, comprises a rotor, defining 
20 rotor poles, a stator defining stator poles, and a set of windings arranged in relation to 
the stator poles to define one or more phases. In a reluctance machine, energization of 
one or more phase windings sets up a magnetic flux in the associated stator poles, 
urging the rotor into a position of minimum reluctance. Timing the sequential 
energization of the windings according to rotor position induces rotor movement. 
25 Switched reluctance machines are well known. More detail is provided in the paper 
'The Characteristics, Design and Applications of Switched Reluctance Motors and 
Drives' by Stephenson and Blake, presented at the PCIM '93 Conference and 
Exhibition at Nurnberg, Germany, June 21-24, 1993 which is incorporated herein by 
reference. As is well known in the art, these machines can be operated as motors or 
30 generators simply by altering the timing of the application of the excitation to the 
phase windings. 

1 
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As explained in the above paper, the method of torque production in a switched 
reluctance machine is quite different from that in conventional machines, e.g. 
induction or synchronous machines, which are operated by rotating waves of magneto- 
5 motive force (mmf) and in which the torque is produced by the interaction of a 

magnetic field with a current flowing in a conductor. Such machines are known as 
'electromagnetic' machines and encompass, e.g., so-called brushless DC machines in 
which the current is in stator coils and the field is produced by permanent magnets on 
the rotor. These machines require the use of permanent magnets. 

10 

By contrast, switched reluctance machines are purely 'magnetic' machines, where the 
torque is produced solely by the magnetic field as the reluctance of the magnetic 
circuit changes. It follows that the methods of controlling the two types of machine 
are quite different, since the control is related to the method of torque production. In 
1 5 general, the control methods used for conventional, sinusoidally fed machines are 
quite inappropriate for switched reluctance machines. 

Figure 1 shows a typical switched reluctance machine in cross section. In this 
example, the stator 10 has six stator poles 12, and the rotor 14 has four rotor poles 16. 

20 Each stator pole carries a coil 18. The coils on diametrically opposite poles are 
connected in series to provide three phase windings. Only one phase winding is 
shown, for clarity. The control of the switched reluctance machine can be achieved in 
a variety of ways. The machine could be controlled in an open-loop fashion, i.e. as 
commonly used for stepping motors. In this regime, the phase windings in the 

25 machine are sent pulses in turn and it is assumed that the rotor lines up with each pair 
of stator poles in turn, i.e. the position of minimum reluctance for that phase which is 
excited. Of course, because the system is open-loop, there are no means of knowing if 
the rotor has moved or not. To remove this uncertainty, it is conventional to use a 
rotor position detection scheme of some sort which provides a signal representative of 

30 rotor position. The excitation can then be applied as a function of the position. Such 
machines are often referred to as "rotor position switched machines". 
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Since current in the windings is relatively easy to measure, closed-loop current control 
is commonly accomplished by monitoring and controlling the energizing current in the 
windings. However, the desired output of the machine is usually torque, position or 
speed, and current has a highly non-linear relationship to all of these. The result is 
that current control techniques generally have inaccuracies in the output, such as 
torque ripple, position error or speed error. 

A typical switched reluctance drive is shown in Figure 2. In this example, the 
machine 36 corresponds to that shown in Figure 1 . The three phase windings A, B 
and C are switched onto a d.c. supply V by a set of power electronic switches 48. The 
moments at which the switches operate are determined by the controller 38, which 
may be implemented either in hardware or in the software of a microcontroller or 
digital signal processor. The firing signals are sent to the switches via a data bus 46. 
Closed loop current feedback is provided by sensing the phase currents by one or more 
current sensors 44 and feeding back signals proportional to phase current. The 
control algorithms often include a proportional (P), proportional-plus-integral (P+I), 
time optimal, feedback linearized, proportional/integral/derivative (PID) function, or 
one of many others as is well understood in the art. It is also common for an outer 
control loop of position or speed to be provided by feeding back a rotor position signal 
from a position detector 40. 

In operation, a current demand i D on line 42 is provided to the controller and this 
regulates the current in the windings, according to the particular control scheme 
adopted, to produce the desired output from the machine. Those skilled in the art will 
be familiar with the many variations of current controllers which exist, each of which 
has its own merits, but all of them suffer from the problems of non-linearity between 
the controlled variable and the machine output described above. 

It has been recognized by the inventor that the more fundamental control variable in a 
switched reluctance machine is the flux which is set up in the magnetic circuit in the 

3 
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machine when a phase winding is energized. The flux is directly responsible for the 
force which acts on the rotor to urge it to a position of minimum reluctance, i.e. to pull 
the rotor round, with respect to the energized stator poles. Embodiments of this 
invention use closed-loop real-time determination and control of flux to achieve much 
better performance from the machine than has hitherto been possible with closed-loop 
control of current. 

In practice, it has been realized, current control is particularly sensitive to variations in 
the machine. In contrast it has been found that flux control is substantially dependent 
only on the flux-causing voltage across the phase winding(s). Thus, it is found to be 
more tolerant of these variations in the machine. Many of the practical uncertainties 
arising from the manufacturing variations to which current control is prone are simply 
avoided. 

In the paper 'Torque Control of Switched Reluctance Drives' by P.G. Barrass and 
B.C. Mecrow, ICEM'96 Proceedings, International Conference on Electrical 
Machines, 10-12 September 1996, Vigo, Spain, Vol 1, pp 254 - 259, incorporated 
herein by reference, there is a proposal to provide a torque control by reference to flux 
linkage reference waveforms using a look-up table that stores fixed values of flux 
ramps for co-ordinates of supply voltage, phase current and rotor position. The flux 
values and co-ordinates are specific to a particular motor. At any instant the pre- 
stored values of flux and torque are chosen from measurements of phase current and 
the stored machine data. There is a fixed relationship between the monitored variables 
and the values of the flux waveforms in the look-up table that are used to produce an 
output for a given motor. 

SUMMARY OF THE INVENTION 

It has been recognized by the inventor that the machine operating characteristic that is 
relatively easily determined, and that relates directly to the flux, is the voltage across 
the phase winding, according to embodiments of the invention. 

4 
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According to embodiments of the present invention there is provided a brushless 
electrical machine comprising: a rotor; a stator; at least one phase winding arranged to 
establish flux in a magnetic circuit in the machine; and means arranged in relation to 
5 the magnetic circuit which are operable to produce a signal indicative of the flux- 
causing voltage across the at least one phase winding. 

According to embodiments of the invention the machine, which can be run as a motor 
or a generator, derives the feedback signal indicative of the flux from the voltage 
10 across the phase winding. In one form the voltage is determined by a transducer 

arranged to measure it directly. A search coil could be used in relation to the phase 
winding, or each phase respectively, which coil will produce a voltage output 
according to the voltage across the phase winding. 

1 5 Also according to embodiments of the invention there is provided a brushless 

electrical drive system comprising: a brushless electrical machine having a rotor, a 
stator and at least one phase winding arranged to establish flux in a magnetic circuit in 
the machine; means for determining the flux-causing voltage across the or each phase 
winding and producing a feedback signal representing the flux-causing voltage; and a 

20 flux controller having an input signal representing the demanded output of the 

machine, which controller is responsive to the input signal and the feedback signal to 
produce control signals for actuating the switch means to control the flux in the at least 
one phase winding. 

25 In a further form, the voltage across the phase winding can be estimated by a state 
observer arrangement, which estimates the flux -causing voltage. 

The voltage across a winding should be integrated to derive the actual flux signal as 
will be explained. A disadvantage of this is that the time between integrator resets 
30 increases with decreasing speed of the rotor, leading to inaccuracies which, of course, 
become intolerable particularly at standstill. Thus, it is also possible to estimate the 

5 
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flux from the current in the phase winding using a current model of the machine, for 
example, according to a state observer arrangement. While there is no integral term in 
deriving the flux from a current model, and it is therefore unaffected at low speeds, it 
makes reference to the angular position of the rotor with respect to the stator and is 
5 therefore less accurate at higher speeds. 

The current model is not sensitive to phase resistance, so the current model can be 
usefully arranged to dominate at "low speed" when little voltage is necessary to 
generate the required flux linkage. If the system is adaptive, the phase resistance can 

10 be more easily estimated at low speed while the observer is relying on the current 
model. At high speed the current model will be less accurate due to angle 
measurement errors being magnified as the speed increases. More applied phase 
voltage is required, according to embodiments of the invention, to generate the 
necessary flux linkage over the shorter phase cycle time. The voltage model is not 

1 5 sensitive to angle measurement error and is less sensitive to phase resistance at high 
speed. The voltage model (the integral with respect to time of the difference between 
the voltage applied to a phase and the voltage drop across that phase resistance) will 
also model complex effects like mutual inductance between phases in a simple 
manner. The total phase flux linkage is derived no matter which phase mmf is 

20 responsible. 

Embodiments of the invention enable direct flux control of the machine which has 
been found to be more accurate and is more amenable to on-line adaptation than the 
current-based control previously used. It uses a real-time determination of flux. 

25 

Preferably, the flux controller further includes means for timing the control signals for 
actuating the excitation means. 

Embodiments of the invention also extend to a method of controlling a brushless 
30 electrical machine having a rotor, a stator and at least one phase winding, the method 
comprising: producing a feedback signal including a part indicative of the flux- 

6 
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causing voltage across the or each phase winding; producing an input signal 
representing the demanded output of the machine; and controlling energization of the 
at least one phase winding in response to the input signal and the feedback signal. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention can be put into practice in various ways, some of which will 
now be described by way of example with reference to the accompanying drawings in 
which: 

10 Figure 1 is a schematic diagram of a known switched reluctance machine; 

Figure 2 is a schematic diagram of a known closed loop current controller for a 
switched reluctance machine; 

Figure 3 is a schematic block diagram of a control system for a switched 
reluctance machine incorporating an embodiment of the invention; 
1 5 Figure 4 is a radial cross-section of another embodiment of a reluctance 

machine according to the invention; 

Figure 5 is a schematic block diagram of a further embodiment of a control 
system in accordance with the invention; 

Figure 6 is a schematic block diagram of part of the system of Figure 5 
20 according to one particular form of the invention; 

Figure 7 is a more detailed schematic block diagram of Figure 6 for a three 
phase machine; 

Figure 8 is a more detailed schematic block diagram of an alternative form of 
Figure 6; and 

25 Figure 9 is a plot of the gain of the current estimator control law with motor 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 



30 



Referring to Figure 3, a first embodiment of a switched reluctance drive system 
according to the invention includes a switched reluctance machine 36 having a 
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laminated stator 10 and a laminated rotor 14 in relation to which is arranged a rotor 
position transducer (RPT) 40. A flux controller 50 produces firing signals for 
actuating power electronic switches 48 for controlling the energization of the phase 
windings A, B and C associated with the stator 10. The schematically illustrated 
5 connection of only one phase winding is shown for the sake of clarity, but each phase 
is independently excitable by actuation of the switches. In an alternative embodiment, 
the excitation for the windings is provided by an amplifier of known type. It will also 
be appreciated that each phase winding will have associated with it a current 
transducer 44. Only one is shown for the sake of clarity. The input 42' to the flux 
10 controller 50 is a signal representing the desired flux v|/ D in the machine to produce a 
desired output. In operating the machine as a motor, the desired output is torque, 
speed or position. In operating the machine as a generator, the desired output is 
electrical power. The flux demand signal can simply be in the form of a scaled 
voltage input representing (e.g.) the torque or speed range of the machine. 

15 

The flux controller of Figure 3 is quite different in operation from the controller of 
Figure 2. Instead of a current demand controlling the output of the machine, the 
output is controlled directly by a flux demand: instead of current being the principal 
parameter fed back to correct the output, flux is now the principal controlled 
20 parameter. This is quite different from all known practicable controllers for switched 
reluctance machines. 



In the flux controller 50 of Figure 3, the flux feedback signal \\r is compared with the 
demanded flux 42' to produce an error signal. The flux controller can follow any one 

25 of a range of control laws, such as proportional, proportional-plus-integral, 

proportional/integral/derivative, time optimal, feedback linearized action, etc. as 
described previously. The machine 36 also has a conventional RPT 40, whose output 
9 provides timing information to the flux controller 50 in conventional manner. 
Alternative embodiments use so-called sensorless position detection systems which 

30 estimate position from other parameters of the machine, as is known in the art. The 
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flux controller 50 may also receive a current signal i from the current transducer 44, 
though this is not used in the manner of a conventional current controller: rather, it is 
used simply to provide overall monitoring of current, to ensure that it stays within a 
predetermined safe level, and/or feed forward I n *R n compensation. The flux 
5 controller 50 produces firing signals for the switching devices 48 that control 
energization of the phase windings. 

An alternative method of deriving the flux signal \|/ is shown in Figure 4, in which 
stator 10' has one or more stator poles of each phase winding 18 that carries an 

1 0 additional winding which acts as a search coil 62. The search coil 62 may be wound 
independently of the phase winding, as shown schematically in Figure 5, or it may be 
wound with it in bifilar fashion. In either case, it will have a voltage induced in it 
when the phase winding is excited which, when integrated with respect to time by an 
integrator 64, gives a signal directly proportional to the flux linking the search coil 62. 

1 5 The integrated signal is fed back to the flux controller 50 as the flux feedback signal 

V|/. 

Figure 4 shows one search coil 62 for the phase winding, but it will be recognized that 
a system with one coil on each pole could be used or, alternatively, the coil(s) could be 
20 placed around some other part of the magnetic circuit, e.g. around the yoke of the 

stator core, on a section between two stator poles, where it is able to sense the flux in 
the magnetic circuit of the machine. 

It will be realized that the embodiments shown in Figures 3 and 4 are illustrative only 
25 and that other methods of producing the flux feedback signal i|/ could be used. The 
main requirement according to embodiments of the invention is that the signal \j/ used 
by the flux controller is proportional to the flux in the magnetic path associated with 
the phase winding being excited. 
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From the above, it will be apparent to those skilled in the art that real-time, flux-based, 
closed-loop control can be applied with benefit to any switched reluctance machine 
with any number of poles and phases. Accordingly, the above description of 
embodiments is made by way of example and not for the purposes of limitation. 

Referring to Figure 5, in another form of the invention a switched reluctance drive 
comprises a switched reluctance machine 1 10 having a stator 112 and a rotor 1 14 in 
relation to which is arranged a rotor position transducer (RPT) 116. The RPT 1 16 
provides rotor position signals G to a flux controller 1 1 8 which produces firing signals 
for actuating a switching circuit 120 for controlling the energization of the windings 
122 associated with the stator. The flux controller 1 18 is also provided with current 
and supply voltage signals i and v from the machine. 

The flux controller must be able to estimate, or 'model' the flux in the machine, given 
the quantities of voltage, and/or current and rotor position. It is possible to achieve 
this in a variety of ways, two of which will be described in detail. Firstly, the flux can 
be modelled from a knowledge of the voltage applied to the phase and the time for 
which it is applied. This is called hereafter the 'voltage model'. It can be refined by 
introducing secondary terms involving current and winding parameters, but essentially 
is driven by the applied voltage. Secondly, the flux can be estimated from a 
knowledge of the current in the winding and the rotor position. This is called hereafter 
the 'current model'. While both of these models can take different analytical or 
numerical forms, preferred implementations will be described below. 

According to this embodiment of the invention, the flux controller 118 uses a state 
observer to model flux linkage, the observer including both a current model and a 
voltage model of the flux. It is found that the current model becomes less accurate at 
high speeds because of increased rotor position sensor error at high speed. The 
voltage model involves cyclical integration of voltage with respect to time and is, 
therefore, susceptible to errors at zero and low speeds at which the reset of the 
integration term occurs increasingly infrequently, leading to drift. Thus, the frequency 

10 
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response of the current model part of the flux estimating system of the embodiment is 
arranged so that it is dominant at zero and low speeds while the voltage model part is 
dominant at higher speeds. 

5 The flux y in a winding is proportional to the integral of the emf, E, applied to the 
winding with respect to time: 



Therefore to control accurately the flux in a switched reluctance machine in real time 
it has been found that it is beneficial to derive a value for the voltage E producing the 
flux \\). To do this, a search coil or bifilar winding can be arranged in relation to each 
phase winding, as described above in relation to Figure 4. 

An alternative approach is to develop an estimate for the voltage E which can be based 
on V, i and a thermal model of, or electrical data for, the winding resistance R. This is 
less intrusive than a search coil or bifilar winding. Since V is the voltage applied to 
the phase winding, i is the phase current and R is the phase resistance, 



20 Although if the resistance of winding R drifts then the flux value drifts, this is still a 

useful method of estimating the flux and, in any case, a thermal model can be set up to 
track any changes in R. However, in many applications it will be acceptable to choose 
a single value of R associated, for example, with a set of standard operating conditions 
or a single operating point. This voltage model will also observe complex effects like 

25 mutual inductance between phases in a simple manner. The total phase flux linkage is 
derived no matter which phase mmf is responsible. 

As with any measurement system involving an integrator, care must be taken either to 
prevent long term drift in the integrator or to avoid depending on its output when drift 




E = V - iR 



so 




[V - iR] dt 
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is likely to be a problem. In this particular case, it is of most concern at or near zero 
machine speed because the integrator has a significant period between integrator resets 
during which drift may occur, and no reset at all at zero speed. 

5 It is therefore preferable not to rely on the estimate produced by the voltage model at 
low speed, but to arrange for a more reliable current model to be used. 

As i|/= J V-iRdt 
10 Flux O = (Ni)/5R (where is the reluctance of the magnetic circuit and N 



is the number of turns in the phase winding) 



d/dt (Flux linkage \|/) = V - iR 



Thus 



15 




and 



i = 9?*<D/AT 




20 



6 = (V/N)-R* ct> *WN 2 , the derivative of the flux estimate with respect to time. 
9?*R varies with rotor angle and temperature 



2 



N 



V is the true control input to the motor. 



25 



N 
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Using a combination of the voltage and current models, a flux estimator equation can 
be written: 



10 



6 = (V/N)-R* 6 *WN 2 +Kobs((x))*(i- O *SRIN) 

Where Kobs is the "observer gain" for the current model which is a function of rotor 
angular velocity co. 

R* 0> *SRIN 2 disappears with the use of a search coil or a bifilar winding. 

Now, Kobs is a time domain non-linear function which is difficult to derive, but at 
high speed the value of Kobs will go to zero, whereas the 6 term contributes less at 
low speed. Kobs becomes large at low speed and dominates the observer as Kobs is 
based in the rotor angle domain. Thus, this embodiment of the invention decouples 
the time and angle domain terms by splitting the estimator between current and 
voltage models and allows each to dominate at different machine speed ranges. 



This basic technique of using both models in a closed-loop controller is shown in 
Figure 6, in which a current model 130 has inputs of phase current i and rotor position 

20 9. The totality of the estimator is indicated in Figure 6. The output of the current 

model 130 is an estimate of flux y 1? which is fed to a summing junction 32 in which it 
is combined with the flux estimate y> to produce a difference value e\\i l which is 
provided to a gain matrix 134. In this embodiment the gain matrix 134 is a 
proportional-plus-integral (P+I) observer gain matrix (error convergence law). The 

25 observer gain matrix (error convergence law) can take other forms, such as a PD, PID, 
hysteretic or predictive controller as will be apparent to the skilled person. 

Block 138 receives inputs of phase voltage V, phase current i and estimated winding 
resistance R and computes the estimate of emf E. (Note that, strictly speaking, this 
30 block is only part of the voltage model). Since E = d/dt (y), this value of E can be 

13 
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summed with the value of d/dt(>/) output from the gain matrix 134, in the summing 
junction 35 to produce d/dt ( y> ). This signal is then supplied to an integrator 136 to 
produce the flux estimate y/ . The flux estimate y> is then applied to a comparator 
140, together with the flux demand signal y/ re f, to produce a flux error signal y/ e which 
5 is applied to control laws 121 to control the switching circuit 120 of Figure 5. 

In Figure 6 the voltage model 138 for the flux-creating emf E is derived from the 
phase voltage, phase current and phase resistance estimate. A thermal model to 
determine the phase resistance R by estimation, or an electrical measurement of phase 
1 0 resistance R can be used or a combination of the two. A convenient form of electrical 
measurement would be by means of a temperature transducer, such as a thermocouple. 
Such a transducer will usually need a filter to attenuate noise that would otherwise 
lead to errors in the measurements. 

1 5 It will be understood by those skilled in the art that Figure 6 represents a controller for 
one phase of the machine. If the machine is polyphase, then each phase will 
essentially have a controller as shown in Figure 6, though in practice there will be 
some sharing of circuitry for economy. 

20 The closed loop current model-based estimator fixes drift at low and zero speeds, 
avoiding distortion of the flux estimate by using the current model in place of the 
voltage model. The P+I gain matrix 134 has a low enough gain and a response 
characteristic such that the flux estimate based on the current model from the matrix 
134 is not influential at high speeds where the voltage model 138 is dominant, whereas 

25 the estimator control law function allows the current model to dominate at low speed. 
Thus, one aspect of the control regime is that the flux estimator hands over between 
the current model and voltage model control techniques according to machine speed. 

The closed loop flux controller will, therefore, have a high "current model gain" at 
30 low speed so the current model dominates and a low or zero current model gain at 

14 
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high speed where the voltage model is more suited (if the current model gain is zero, 
the integrators should be reset at known "zero flux" intervals). The current model (130 
of Figure 6) could be a look up table or a functional description in the form of an 
algorithm. The advantage of the functional description is that it is amenable to 
adaptation via an adaptive mechanism. The functional description can be linear in the 
coefficients that change via adaptation. Even without the adaptive mechanism, the 
system uses both the current and voltage models to estimate flux linkage. If the 
system is adaptive, at low speed the phase resistance will be corrected, whereas at high 
speed the "current model" and potentially phase resistance will be corrected. The 
result is essentially no interaction of the adaptive mechanism with the flux observer. 

To summarize: the voltage model becomes less and less accurate as motor speed 
decreases toward zero because: 

• The open loop integration process will experience drift and random walk 
(uncertainty growing with time as a result of integrating a signal containing noise) 
for longer intervals between resets. 

• At low speed the applied voltage is dropped primarily across the phase resistance 
which changes with temperature. The effects of the angle varying inductance will 
be very small at low speeds resulting in a poor signal to noise ratio. 

• If a low pass filter is used to approximate an integrator at high speeds, the error 
between low pass filtering and true integration increases as the motor speed 
decreases. For zero mean a.c. signals a low pass filter can be used alone. If the 
signals have a d.c. content there must be a reset function or k obs must have a 
sufficiently high gain. 



At high speed the voltage model is accurate and at the same time very simple. The 
voltage model is more accurate at higher speeds because: 

• The time between resets is much shorter resulting in less drift and random walk. 

• At high speeds the voltage is dropped primarily across the angle varying 

30 inductance which causes less estimation error due to phase resistance changes. 
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• In the case of using approximate integration (e.g. a low pass filter) the difference 
between true integration and low pass filtering becomes negligible at high speed 
for a.c. signals. 

5 The voltage model is a computation of V- i*R (the flux-causing emf E) from which 

the flux is derived by reset integration or low pass filtering. All of the complex mutual 
effects are included in this calculation no matter which phase mmf is driving the flux. 
A complex look-up table, or calculation of a complex function involving all motor 
phase currents and rotor position, can be very computationally intensive requiring ever 

1 0 more computing power as the motor speed increases. This is a further reason why the 
current model may not be used at high speeds. Thus, the voltage model can be used to 
perform the flux estimate in isolation at high speeds. The voltage model does not use 
the shaft position to instantaneously estimate flux linkage. The increasingly inaccurate 
shaft position measurements/estimates with increasing motor speed (calculation 

1 5 delays, signal filtering delays, loss of resolution, etc.) are not of concern to the voltage 
model. 

In situations in which less accurate control of the switched reluctance machine is 
acceptable, or in which flux control is only required at sufficiently high speeds, it is 
20 possible to use only the voltage model and not to apply the current model in 

controlling flux . However, this does not provide the accuracy at zero and low speeds. 

Figure 7 is a more detailed three phase example of the flux control system described in 
Figure 6. Like reference numerals have been used for equivalent parts in Figures 6 
25 and 7. It will be noted that each phase has separate inputs and outputs in respect of the 
various values of current, voltage and flux. 

The three phase SR example of the embodiment described by Figure 7 includes a 
closed loop flux observer having the voltage model and the current model 130 of the 
30 SR machine phase flux, which predict the flux for each phase as a function of rotor 
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position 0. The respective models are summed into the observer at summing junctions 
135', 135" and 135"'. 

As stated above, the current model can take many forms, from a look-up table to a 
functional description of the current in the form of the real-time processed current 
model. Below are examples of current models that can be used: 

Current Model #1 can be presented in the form of a current-product relationship: 

¥in= (aii n + a2i n 2 + a 3 i n 3+ ) 

+ (b!i n + b 2 i n 2 + b 3 i n 3 + )cos(e en ) 

+ ( C1 i n + c 2 i n 2 + c 3 i n 3 + )sin(9 en ) 

+ (d 1 i n + d 2 i n 2 + d 3 i n 3+ )cos(29en) 

etc. 

Where 9 is the electrical or mechanical angle and n is the phase number. 0 < 9 en < 
2n for phase n. The mechanical angle is used in the model if electrical sub-harmonics 
are to be included in the model. At any angle, flux is a function of current. However, 
the coefficients of such a relationship are dependent upon angle. Such a relationship 
can be described using the current-product relationship presented above. 
Mathematically such a relationship is, in the limit, true to arbitrary accuracy due to the 
nature of the trigonometric and polynomial functions used by virtue of their ability to 
approximate any continuous function. Electrical sub-harmonics are typically the result 
of manufacturing variations in a SR machine. The form of Current Model #1 allows 
each harmonic of a phase flux linkage to change in an arbitrary "smooth" non-linear 
fashion as a function of that particular phase current. Mutual effects are not 
considered in Current Model #1. 

Current Model #2 uses a hyperbolic tangent form. The generation of flux from current 
occurs through the physical medium of the motor iron. This material exhibits a 
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relationship between current and flux which can be closely approximated with the 
hyperbolic tangent function. Such an expansion is also well behaved with respect to 
out of range values: it exhibits a limiting behavior. Hence, out of range values for the 
object variables do not produce abnormal or pathological behavior in the observer 
5 model. 

The hyperbolic tangent function may take the form: 

Vin = (Vn + a l tanh (in) + a 2 tanh ( { n 2 ) + • • • + a p tanh ( [ n v )) 
+ (bi tanh (i n ) + b 2 tanh (i n 2 ) + . . . + b p tanh (i n P)) * cos (6 en ) 
10 + (c i tanh (i n ) + c 2 tanh (i n 2 ) + . . . + c p tanh (i n P)) * sin (0 en ) 

+ (di tanh (i n ) + d 2 tanh (i n 2 ) + . . .+ d p tanh (i n P)) * cos (29 en ) 
+ (ei tanh (i„) + e 2 tanh (i n 2 ) + . . . + ^ tanh (i n p )) * sin (20 e „) 
+ ... 
alternatively, 

15 \|/ in = (a 0 i n + a! tanh (i n ) + a 2 tanh (2*i n ) + . . . + a p tanh (p*i n )) 

+ (bj tanh (i n ) + b 2 tanh (2*i n ) + . . . + b p tanh (p*in)) * cos (0 en ) 
+ (ci tanh (i^ + c 2 tanh (2*i n ) + . . . + c p tanh (p*i n )) * sin (0 en ) 
+ (di tanh (i n ) + d 2 tanh (2*i n ) + . . .+ d p tanh (p*i n )) * cos (29 en ) 
+ (eO tanh (i n ) + e 2 tanh (2*i n ) + .. . + ep tanh (p*i„)) * cos (2e en ) 

20 +... 

or any mixture of the above two model forms. The above two model forms do not 
consider mutual effects. 

Conveniently, the complicating factor, the effect of mutual inductance between 
25 phases, can be addressed by the use of look-up tables containing measured data or 
another functional form. 



18 



Attorney Docket No.: K315.106.101 

Returning to the K G b s observer gain matrix 134 of Figure 6 or 7 it can be any function 
that ensures that during low and zero speed operation goes to zero as time (t) goes 
to infinity (ei|/ n is the estimation error between the current model estimate of phase 
flux linkage, \jr i n , and the complete estimate of phase flux linkage, y/ a ). During high 
5 speed operation K c b s may be set to zero provided the integrators are reset. The low 
pass filter option may be used with zero mean AC signals for the flux-causing voltage. 

The example for a 3-phase machine shown in Figure 7 is based on the diagonal matrix: 



~K4a>) 0 0 
0 Ki(oS) 0 
0 0 Ks(co) 



where K\ , K2 and K3 are dependent on the angular speed (go) of the rotor or are 
constant gains. Constant gains would be used only in a simplified implementation 
with reduced performance. K 0 ^ s can be diagonal in this version of the flux observer 
1 5 current model to simplify calculation of K Q b s . 



Note that mutual effects are considered in the flux observer models of Figures 6 and 7. 
The current model 130 estimates each phase flux linkage as a function of every phase 
current and angle. In addition the total phase (coil) voltage is used by the voltage 

20 model portion of the flux observer. The integrator for the or each phase in the voltage 
model flux observer is reset during known periods of zero (or as near as possible zero) 
flux in the phase cycle. At high speeds each integrator could also be adapted to act as 
a low pass filter (with, for example, time constants about 10 times longer than the 
fundamental frequency of the motor phase switching period) when the supply to the 

25 machine has a zero mean a.c. characteristic or a unipolar signal is transformed to a 
zero mean a.c. representation. 
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Writing current model #1 in summation form: 



harm_max poly_max 

l}/_I_x_cos(i_x,e,harm_max,poly_max) =Z £ L_x_i_cos n m (i_x) m cos(n0) (1) 



n = 0 



l[/_I_x_sin(i_x,8,harm_max,poly_max) 

10 

\|/_x-I_x_tot=l[/_I_x_cos(i_x,e ,harm_max,poly_max)+V|/_I_x_sin(i_x 9,harm_max,poly_ max) (3) 

where harm_max is the highest harmonic of the phase flux linkage profile considered 
1 5 and polyjmax is the highest polynomial power considered. This is the expression for 
current model # 1 shown above, but in summation notation. 

In a symmetric machine which does not exhibit significant remanence properties due 
to effects such as lamination grain orientation in the material of the rotor and the 
20 stator, V|/_I_x_tot should have only odd polynomial powers. In other words the self- 
induced flux linkage should have the same sign as the current that produced it. 

The effect of mutual inductance or coupling between the phases, in addition to the 
phase flux linkage generated by the phase current and its self-inductance, must be 
25 considered when determining the total phase flux linkage. The steel of a switched 
reluctance machine is a non-linear medium for flux. Thus, a general form of model 
for mutual flux linkage between phases is required. To take the example of a three 
phase machine, the general form of expression above can be modified to be written as: 



harm_max poly_max 

= S Z L_x_i_sin n , m (i_x) m sin(n9) (2) 

n=0 m=l 
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y_x_I_x_I_y_I_z_cos(i_x,i^4_z,6,m^ = 



ESSE ^jiJ_xJjyJ_z_cos k . mkn i_x k iy n i_z b cos(nQ) (4) 



\j/_x_I_x_I_y_I_z_sin(i_x4_y4_z,B,m_harm_max,poly_x_max,poly_y_max,po - 

10 poly_x_max poly_y_max poly_max_z m_harm_max 

I Z Z I v_x_i_x_i_y_i_z_sin k>m ,i 1 , n i_x k i_y m i_z h sin(ne) (5) 

If no magnets nor any significant remanence effects are present (which is generally the 
1 5 case for the switched reluctance machine): 

y_x_I_x_I__y_I_z_cos 0 ,o,o,n=0 and v|/_x_I_x_I_y_I_z_sino >0j o,n=0 (6) & (7) 

Thus the total flux linkage in any phase can then be written as: 

20 

v|/_x_tot=\|/_x_I_x_I_y_I_z_cos + vi/_xJ^_x_I_y_I_z_sin (8) 

Since equations (4) through (7) include both the mutually induced flux linkages, which 
may be affected by mmf polarity-dependent magnetic path changes, and the self- 
25 induced flux linkages, no generalization can be made concerning the polynomial 
exponents as can be made for the self-induced flux linkages in isolation. 

Equations (4) and (5) above include both the self-induced flux linkage and the 
mutually induced flux linkage in a given phase. If the self-induced flux linkage is to 
30 be separated from the mutual flux, equations (4) and (5) can be rewritten with different 
constraints as: 
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m_x_I_x_I_y_I_z_cos(i_x,i^,iz,9,mJiarm_max,poly_x_m = 
Z Z Z E m_x_i_x_i_y_i_z_cos kjm , h , n i_x lc i_y m i_z h cos(n.9) (9) 

k=0 m=0 h=0 n=0 

m_X_I_X_I_y_I_Z_COS 0 p0 l y x max,0 ,0 n = 0 ( 1 0) 

m_x_I_x_I_y_I_z_sin(i_x4_y,i_z,8,m_harm_max,poly_x_max,poly_y_max,poly_max_z) 



poly_x_max poly_y_max poly_max_z m_harm_max 

Z Z Z Z m_x_i_x_i_y_i_z_cos kjmihin i_x k i_y m i_z h sin(n.9) (11) 



m_X_I_X_I_y_I_Z_COS 0 po ly_x_max,0 ,0,n = 0 

The total mutually induced flux linkage can then be written as: 

20 m _ x _J_ x _I_y_I_ x tot = m_x_I_x_I_y_I_z_cos + m_x_I_x_I_y_I_z__sin (12) 

Equations (9) through (12) include all of the flux linkage generated in phase x that is 
not due to I_x alone, in other words, all of the mutual inductances. Since equations (9) 
through (12) describe the mutually induced flux linkages which may be affected by 
25 mmf polarity-dependent magnetic path changes, no generalization can be made 
concerning the polynomial exponents as can be made for the self-induced flux 
linkages in isolation. 

The total flux linkage in a given phase would then be written as: 
30 V_ x _tot(i_x_i_y_i_z,9) = v)/_x_I_tot+m_x_I_x_I_y_I_z_tot (13) 



Using equations of the form of (4) through (7) the current dependence of the flux 
linkages can be generalized as: 
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\|;_xj_xj_y_l_z_cos(i_x4_y4_z,0,m_harm_max,poly_x_max,poly_y_max,poly_max_z) = 
X Z Z Z V_x_i_x_i_y_i_z_cos k , m , h , n f(i_x k ).f(i_y m ).f(i_z h )cos(n.e) 

k=0 m=0 h=0 n=0 (U) 



V(/_x_I_x_I_y_I_z_sin(i_x4_y4_z,©,ni_harm_max,poly_x_max,poly_y_max,poly_max_z)= 
Z Z Z Z V _x_Lx_i_y_i_z_sin k , nl , h , n f(i_x k ).f(i_y m ).f(Lz h )cos(n.e) 

k-0 m=0 h=0 n=0 (15) 

If no magnets nor any significant remanence effects are present (which is generally the 
case for the switched reluctance machine): 



20 



^_x_l_l_yj_z_cos 0?0)0n =0 (16) 

Y_ x JLLyJL z _ sil H),0,0,n =0 ( 17 ) 

One example of the form of the functions of current seen in equations (14) and (15) is 
the hyperbolic tangent: 



/(i_x k ) = tanh 
/0_y*) = tanh 
/(i_z k ) = tanh 



f ?=£ Y 

^ norm _ i _ x _ factor ) 



norm _ i _ x _ 


factor ^ 






norm i y 


factor 


i z 





norm _i _z _ factor 



(18) 
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The hyperbolic tangent is just one of a large number of functions that could be chosen. 
Others include spline descriptions, exponential, etc. The important property of the 
equations described is that the equations are linear in the coefficients for each poly- 
harmonic product term (each poly-harmonic term in equations (14) and (15) are 
5 linearly dependent on the appropriate H/_x_l_x_l_y_l_z_sin_k_m_h coefficient). The 
equations being linear in their coefficients allows for least square parameter 
estimation. Least squares-based parameter estimation is a well proven method for 
adapting system models in an embedded application. 

10 The creation of functional descriptions of the system which are linear in coefficients 
allows for the application of the theory of linear systems, i.e. parameter estimation 
techniques, which are well understood and reliable. 

For the linear magnetic case, equations (14) through (17) reduce to 

il/_x_I_x_I_y_I_z_cos(i_x4_y4_z8,m_hami_max,poly_x_max,poly_y_max,poly_max_z) = 
2-. 2-> 2- Z, A-_x_i_x_i_y_i_z_cos k h>n (i_x k ).0_y m )-(Lz b )cos(n.9) 

k=0 m=0 h=0 nK, (21) 

\j/x_I_x_I_y_I_z_cos(i_x,i_y,i_z9,m_harm_max,poly_x_max,poly_y_max,poly_max_z) = 

1 1 1 m_harm_max 

IEEE VK_x_i_ X _i_y_i_z_sin k , m , h , n (i_x k ).(i_y m ).(i_Z h )si n (n.0) 
25 M M »-° (22) 

If no magnets nor any significant remanence effects are present (generally the case for 
the switched reluctance machine) 

30 
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y/ _x I _x I __y _I _z _ cos o.o.o. « 

Y _X _I _X _1 _y _I _Z _ COS 0,1,1. n 

y/_x_I_x_I_y_I_z_ cos 1,0,1, n 
_x _I _x _I _y _I _z _ cos 1,1,0, n 
y/ _x _I _x _I _y _I _ z _ cos 1,1,1, n 



I _x _I _y _I _z _sino, 
I _x I _y _I _z _sino, 
_ / _x I _y I _z _ sin i, 



smi, 



Or in simplified form equations (21) through (24) become: 



10 



>i_x_I_x_I_y_I_z_cos(i_x,i_y4_z,9,harm_max_xx,harm_max_xy,harm_max_xz) = 



Yi A._x_i_x n cos(n.9) + i_y X A__x_i_y n cos(n.9) + i_z X \_x_i_z„cos(n.9) (25) 



A._x_I_x_I y_I_z_sin(i_x,i_y4_z ) 9,harm_max_xx,harm_max_xy,harm_max_xz) = 

15 haml „ m ax.x X hann^x, ha„x^ 

i_x S ^_x_i_x n sin(n.6) + i_y 2 \_x_i_y n sin(n.9) + i_z 2 A,_x_i_z n sin(n.8) (26) 



The total flux linkage can be written as: 

20 

A_x_tot=^_x_I_x_I_y_I_z_cos+^_x_I_x_I_y_I_z_sin 



In both the linear and non-linear examples for a three phase SR machine, model forms 
have been identified that describe the relationship between all three phase currents, 
25 
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angle and phase flux linkage. The expressions described in equations (4) through (27) 
represent current model forms that describe a three phase SR machine, including 
mutual effects. The expressions described above in equations (1) through (3) maybe 
used as the current model 130 of Figure 6 and equations (4) through (27) may be used 
5 as the current model 130 of Figure 7. The model forms described by equations (1) 
through (27) are linear in the model coefficients and have been verified in practice. 
Being linear in the model coefficients allows the model forms described in equations 
(1) through (27) to be tuned via a Recursive Least Squares-based adaptive mechanism. 

10 Figure 8 illustrates an alternative flux observer. In this embodiment a matrix K 0 234 
is a function of phase currents ij, i^, 13, mechanical rotor angle 0 and rotor angular 
velocity g> producing the differential flux estimate output. The flux estimate output of 
the matrix K 0 for each of the three phases is summed with the flux-causing voltage, as 

dy> 

before, to derive a value of the rate of change of flux with respect to time ^ f° r 

15 each phase Thus, the integral of this, from integrator 136, will be the flux observer 
value y/ for each phase. 

In this embodiment, an inverse current model 230 describes each phase current, 
( ^ , i 2 , i 3 ) , as a function of each phase flux estimate, (yr j, y7 2, W 3) and the rotor angle, 
20 9. Values of each phase current are summed with the negative of the output of the 
inverse current model 230, ( i { , i 2 , i 3 ) , in a comparator 232 to produce an error signal 
si n . The output of the inverse current model 230 is derived from the value of the three 
observer fluxes ( y/ \, y/ 2, $3) and rotor angle (9). 

25 This embodiment uses the inverse of the flux-linkage current model that was described 
earlier in conjunction with Figures 6 & 7. The inverse current model is needed (phase 
current estimate as a function of the phase flux or flux linkages and rotor angle) 
because the measured phase currents are compared with the estimated phase currents 
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as shown in Figure 8. In this form the observer is in the classic Luenberger form (see 
"An Introduction to Observers", Luenberger, D.G., IEEE Translations on Automatic 
Control, Vol. AC-16, No. 6, December 1971, pp. 596-602), incorporated herein by 
reference. This flux observer solves the equation: 

5 

ij> =V-i(^e)R + K o [i( V ,G)-i(^,0)] 



where i is calculated from the inverse flux linkage, y/ , ^ , y/, i, V are all nxl vectors 
where n is the number of motor phases. The current model observer gain Kq is based 
10 on a nxn matrix for an n-phase machine. 

From the derivation (not shown) of the observer gain matrix it can be shown that: 

K Q (9,0,^,12,13) = - R n " K(©) / Gi,i 2 ,i3,e) 

15 In matrix form: 







0 


0" 




y,(i\,0) o 


0 


Ko = - 


0 


R„ 


0 


-K((o) 


0 f,(h,6 


) o 




0 


0 


Rn 




0 0 





ignoring mutual inductance. More generally, in the linear magnetic case where mutual 
20 inductance is included: 



Ko = 



Rn 


0 


0 


0 


Rn 


0 


0 


0 


Rn 



-K(co) 



M0) MO) MO) 
MO) f*(6) MO) 
MO) M0) M0)_ 



Two other general, non-linear, cases can be similarly defined: 

25 
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'Rn 


0 


0" 






+ fi^{u,i 




Ko= - 


0 


Rn 


0 


-K(<D) 


MM) 








0 


0 


Rn 






+ fi 3a (u,l 





'Rn 


0 


0" 


0 


Rn 


0 


0 


0 


Rn 



-K(co) 



„ (h,i2,h, 
¥«- — 



, (l../2.»,0) 



, (h,i2,i3,0) 



In these expressions mutual inductance is included. It will be appreciated that the 
form of the expression depends on how the total flux model is defined. The K(co) 
function has a frequency response that is typically similar to that in Figure 9, such that 
the K(co) function dominates at zero and low speeds, but the voltage model dominates 
with increasing speed. Typically, the response of the current estimator will begin to 
tail off at about 10% to 20% of base speed, i.e. the highest speed of maximum torque 
output of the machine. 

The invention provides a flux control technique for electronically switched brushless 
machines of various types, but particularly switched reluctance machines. It produces 
a real-time signal estimate of the flux in the magnetic circuit of the machine to which a 
reference value of flux can be compared to control the machine output. Preferred 
embodiments take account of the differing considerations at high and low speeds by 
using a current model and a voltage model of the machine flux according to machine 
speed. 

It will be apparent to the skilled person that various modifications and changes can be 
made to the specifically disclosed embodiments without departing from the invention. 
The invention is to be limited only by the spirit and scope of the following claims. 
28 



